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ABSTRACT 

Stars form by the gravitational collapse of interstellar gas. The thermodynamic response of the gas 
can be characterized by an effective equation of state. It determines how gas heats up or cools as it 
gets compressed, and hence plays a key role in regulating the process of stellar birth on virtually all 
scales, ranging from individual star clusters up to the galaxy as a whole. We present a systematic 
study of the impact of thermodynamics on gravitational collapse in the context of high-redshift star 
formation, but argue that our findings are also relevant for present-day star formation in molecular 
clouds. 

We consider a polytropic equation of state, P = /cp r , with both sub-isothermal exponents T < 1 
and super-isothermal exponents T > 1. We find significant differences between these two cases. For 
T > 1, pressure gradients slow down the contraction and lead to the formation of a virialized, turbulent 
core. Weak magnetic fields are strongly tangled and efficiently amplified via the small-scale turbulent 
dynamo on timescales corresponding to the eddy-turnover time at the viscous scale. For T < 1, on 
the other hand, pressure support is not sufficient for the formation of such a core. Gravitational 
contraction proceeds much more rapidly and the flow develops very strong shocks, creating a network 
of intersecting sheets and extended filaments. The resulting magnetic field lines are very coherent 
and exhibit a considerable degree of order. Nevertheless, even under these conditions we still find 
exponential growth of the magnetic energy density in the kinematic regime. 



1. INTRODUCTION 

Star formation occurs under a large variety of condi- 
tions, ranging from Milky Way- type molecular clouds, 
starburst galaxies, environments of supermassive black 
holes to the primordial Universe. Environmental con- 
ditions such as the ambient radiation field, the flux of 
cosmic rays, and the overall metallicity of the system, 
determine the thermodynamic response of the gas to 
compression. At low metallicity, for example, the gas 
generally cools less efficiently, increasing the characteris- 
tic temperatures and changing the slope of the effec tive 
equation of state ( Omukai et al.||2QQ5 Omukai 2012]). In 
starburst galaxies and close to supermassive black holes, 
the gas is exposed to UV and/or X-ray radiation fields, 
providing addition al sources of heating , ionization and 
photodissociation ( Malon ey^ et al.||1996| |Kaufman et al.| 
1999; [Meijerink fc Spa ans 2005]). In starburst galaxies, 
the high amount ol cosmic rays may provide a tempera- 
ture plateau at high densities, again having a potentially 
strong impact on the star formation process ( Papadopou- 
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change significantly depending on the conditions of the 
environment. In 3D simulations, the effects of such vary- 
ing equations of state have however only been explored 
in rather specific situations. For instance, strong ef- 
forts were dedicated to the modeling of primordial star 



formation dAbel et al. 20021 


Bromm & Larsonl 2004 


Yoshida et al.||^00«| Clark et aJ 
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mation in metal-poor environments (e.g. Pop eke et al. 
2012] |Glover fc Claxk| |2012b| |Latif et al.| \20V2\ and 



( Klessen et al. 



2007 



orma- 



circum- nuclear starbu rst regions 
Hocuk fc Spaans|2011 ) . For the modelling ol star 
tion in Milky Way-type molecular clouds, different ap- 
proaches exist, rangi ng from piecewise-polytropic equa- 



tions of state (e.g. Bonnell et al. 



2008) to the solution ol the actual heating -cooling bal- 
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In such a variety of different conditions, it is central to 
obtain a more general understanding of how the equa- 
tion of sta te influences stellar birth. For example, |Lar-| 

;ntial 



son 



(2005) provided a detailed discussion of the potential 
role ol the equation of state on the fo rmation of fila- 
ment s. Based on numeri cal simulations, Li et al. (2003) 
and Jappsen et al. (2005 ) studied the impact on the over- 
all fragmentation behavior of the gas, but they have not 
modeled the process of protostellar collapse and stellar 
birth on small scales and neglected magnetic fields. Ana- 
lytic theories suggests a dependence of the l ocal collaps e 
timescale on the equation of state (see, e.g. Yahil| [l983). 

The gas that formed the first stars and p rotogalax- 
ies initially collapses sub-isothermally ( |Scalo fc Biswas 
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2002) and reaches a phase of super-isothermal collapse 
when the gas becomes optically thick. In between these 
two density regimes, the behavior can be very compli- 
cated and chan ge multiple times, depending on the ex- 
act metallicity ( |Omukai et al.|[2010| ). In this letter, we 
systematically explore the impact of different equations 
of state on protostellar collapse and the star formation 
process in the context of primordial star formation using 
3D simulations. We consider weakly magnetized collaps- 
ing minihalos and explore the effect of the equation of 
state on the general morphology, the turbulent proper- 
ties, and the structure and growth of the magnetic field. 
Particular attention is paid to magnetic field amplifica- 
tion processes, which may play an important role even 
at ea rly cosmic times (Schleicher et al. 2010; Sur et al. 
' 2010} iFederrath et al.||2011b[ |Schober et al.||2Ul2fr ~ 

In the following, we outline our numerical approach 
in Section [2j and present our results in Section [3j We 
summarize our findings and discuss their implications in 
Section H 

2. NUMERICAL SIMULATIONS 

We carry out numerical simulations with the adap- 
tive mesh (magn eto-) hydrodynamics code FLASH (|Fryx- 



ell et al. 2000[). We use a positi ve-definite, M USCL- 



Hancock, Kiemann scheme HLL3R (Waag an et al.|20lT ) 



to solve the magnetohydrodynamic equations. Because 
our aim is to study the details of the dynamical evolution 
during protostellar collapse, we adapt ively increase the 
numerical resolution as the density increases and ensure 
that the local Jeans length is always resolved with at least 
64 cells. Previous studies showed that at least 32 grid 
cells per Jeans length are required to resolve vorticity and 
turbulent magnetic- field amplification o n the Jeans scale 
( |Sur et al |2010l jFederrath et al.|2011b| , so we will safely 
capture these effects in our simulations. However, the 
amplification of small-scale magnetic fields and vorticity 
occurs on the viscous scale, where the non-linear term 
in the induction and vortic ity equations lead to an expo- 
nential amplification (e.g., |Federrath et al.||2011a| ). The 
viscous scale in real systems, however, is typically much 
smaller than can be resolved with any current simulation 
technique. Thus, even though the total kinetic energy on 
the Jea ns scale is converged w ith a Jeans resolution of 
32 cells ( Federrath et al.|2011b ), vorticity production on 



smaller scales is not. The amount of vorticity produced in 
simulations like the ones discussed below increases with 



at the same 


time (Sur et al.|2010 Federrath 


Turk et al. 


2012 


). Furthermore, baroclin- 



is absent in our barotropic simulations, can additionally 
contribute to vorticity production. Thus, the magnetic 
field and vorticity produced in our simulations is a lower 
limit of the actual amount produced in real systems. 

For our initial condition we use the central region of 
a primordial minihalo representative of those in which 
the first stars form. The gas in such a minihalo consists 
purely of Hydrogen and Helium, and at the central den- 
sities considered here (10 _15 gcm 
isothermal equation of state (T 
metal-free case. To allow a direct comparison between 
the simulations and a systematic study, we adopt the 
same initial conditions for all equations of state. The 



3 ) follows an almost 
1.1) in the purely 



gas distribution is smooth and centrally concentrated, 
and the gas temperature is high, meaning that the inter- 
nal turbulence is sub- to transsonic. Intrinsic magnetic 
fields in primordial minihalos are thought to be weak, 
but we show here that very small initial magnetic fields 
are exponentially amplified on timescales shorter than 
the collapse. 

Our initial conditio n was extracted fro m the cosmo- 
logical simulations of Greif et al. (2011). We select a 



typical minihalo and extract its central region on scales 
of the Jeans length, corresponding to 8000 AU, with a to- 
tal mass of 101.3 M and a mean temperature of 710 K. 
The initial peak density of the simulation is 2.8 x 10 -15 g 
cm -3 . The original simulation w as carried out using 



the Lagrangian mesh code Arepo (Springel 2010) and 
remapped onto a cartesian grid with a resolution of 256 3 
grid cells to generate the initial conditions for the resim- 
ulations discussed here. 

The initial, turbulent magnetic field has a power-law 
spect rum, Pb(/c) oc k 3 ^ 2 on large scale s ([Kazantsev 
1968) , typical of the small-scale dynamo (Brandenburg 
& Subramanian| [2005 ) and peaks on a scale correspond- 
ing to 20 grid ceil s. The peak scale was set such, because 



Federrath et al. (2011b) found a range of scales corre- 
sponding to 20-30 grid cells, on which the magnetic field 
grows fastest due to small-scale dynamo action. Since we 
want to address potential formation of coherent large- 
scale magnetic fields in this letter, we chose the lower 
limit in order to ensure that we start with a truly small- 
scale field. Below this resolution limit, the spectra drop 
with Pb(/c) oc k~ 4 . The simulation results are, however, 
insensitive to the details of this drop in the dissipation 
range. This initial, turbulent magnetic field serves as a 
sensible initial condition for a field that was previously 
ampli fied by small-scale dynamo action in a collapsing 
halo (IFederrath et al.||2011bl |Xu et al.||2011( iTurk et al. 



2012)! In all simulations we have set the amplitude ol 
the random magnetic field to an initial value of 10 -9 G, 
which is too small to reach the regime where the small- 
scale dynamo saturates during the collapse but allows us 
to study its exponential growth in the kinematic regime 
of the turbulent dynamo. 

To cover the range of possible thermodynamic re- 
sponse behavior of star-forming gas at different metal- 
licit ies and in different environments, we adopt a sim- 
ple polytropic equation of state, P = kp F . We have 
run four collapse calculations with different polytropic 
indices T = 0.7,0.9,1.1,1.2. A value T < 1 corresponds 
to a strong cooling regime as provided, e.g., by molecular 
or atomic line emission at all metallicities or the coupling 
to dust at low metallicites (e.g. |Klessen et al.|2012[ ), while 
T > 1 corresponds to a heating regime as realized when 
the gas becomes optic ally thick or couples to the dust at 
high metallicities (e.g. |Glover fc Clark|[2012a[ ). 

Our primary motivation is the study ol gravitational 
collapse and stellar birth in the early Universe, where 
magnetic fields are initially weak and their subsequent 
amplification is of interest. However, the insights we gain 
from altering the equation of state apply to a range of 
different scenarios and may also be relevant to the dy- 
namics of nearby star-forming clouds in the Milky Way 
at present days. 

3. RESULTS 
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Fig. 1. — Magnetic field and density structure for the four runs as function of time. Rows are different polytropic indices (T = 0.7,0.9, 
from top to bottom), columns are different times (time advances from left to right). The magnetic field vectors have been rescaled for 
plotting by p 2//3 . The field strength mentioned in the images corresponds to an arrow of the length given in the legend. The green circle 
marks the Jeans volume. 
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Fig. 2. — Magnetic field and density structure for the four runs as function of time. Rows are different polytropic indices (T = 1.1, 1.2, 
from top to bottom), columns are different times (time advances from left to right). The magnetic field vectors have been rescaled for 
plotting by p 2 / 3 . The field strength mentioned in the images corresponds to an arrow of the length given in the legend. The green circle 
marks the Jeans volume. 
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The collapse dynamics is remarkably different depend- 
ing on the polytropic exponent. To illustrate this point, 
we show in Figures [T] and [2] density slices and magnetic 
field vectors for the four simulations (rows) at three dif- 
ferent times (columns). The density structure and mag- 
netic field topology changes significantly when the poly- 
tropic exponent shifts from super-isothermal (r > 1) to 
sub-isothermal (r < 1). 

In the super-isothermal case, strong pressure gradients 
are able to build up. They slow down the gravitational 
collapse and the flow tends to virialize in the central core. 
The accretion o nto this core is able to generated strong 



turbulence (e .g. Klessen & Hennebelle 



et al. 2011b). In the kinetic weak-held regime studied 



2010 Federrath 



here, the magnetic field is strongly tangled and efficiently 
amplified via the small-scale dynamo on timescales corre- 
sponding to the eddy-turnover time at the viscous scale 
( |Kazantsevj|1968| | Schober et al. 2012[ ). In contrast to 
that, sub-isothermal collapse exhibits a markedly differ- 
ent dynamical behavior. In the absence of strong (stabi- 
lizing) pressure gradients, gravitational contraction pro- 
ceeds much more rapidly. As a consequence, the flow 
develops very strong shocks, which create a network of 
intersecting sheets and extended filaments rather than 
one central virialized core. The resulting magnetic field 
lines are very coherent and considerably more ordered 
than in the case T > 1. 

The two different regimes can also be seen in the ve- 
locity structure. Figures [3] and [4] show the divergence 
of the velocity field divv and the z-component of the 
vorticity (curl v) 2, both normalized by the local freefall 
time tff, for one snapshot of each simulation. The sim- 
ulations are compared at similar mean densities in the 
Jeans volume. Clearly, the velocities have significantly 
more small-scale structure in the turbulent core of the 
super-isothermal than in the sub-isothermal cases. It is 
interesting to note that even for T < 1 we find both 
compressive and solenoidal turbulent modes close to the 
curved shock structures (see Fig.j3|, but very little veloc- 
ity structure in the infalling regions. This behavior can 
be explained by vorticity production in curved sh ocks 
flMee fc Brandenburg|[2006[ [Federrath et al.||2011aj ) and 
by viscous interactions at density gradients, efficiently 
operating even in regimes of primarily compressive driv- 
ing of turbulence, as during gravitational contraction 
studied here. 

We now turn our attention to the magnetic field am- 
plification. Indeed, the effect of the two different dy- 
namical regimes is clearly visible in Figure [5] We plot 
the turbulent velocity dispersion a Vl the magnetic field 
amplification normalized by the contribution from flux- 
freezing 5/p 2 / 3 , the ratio of the fluctuating parts of the 
kinetic energy and magnetic energy E f mag / E^ in , as well 
as the turbulent Mach number M 1 for the four simula- 
tions. All quantities are averaged over the Jeans volume 
during the collapse and plotted against the mean density 
in the Jeans volume. Contributions from a mean velocity 
or magnetic field have been subtracted to calculate the 
fluctuating components. The growth of the turbulent ve- 
locity dispersion is much weaker for the sub-isothermal 
collapse than in the super-isothermal case. Regarding 
the evolution of the magnetic field, the ratio B/p 2 ^ 3 is in- 
creasing for all four simulations, confirming the presence 



of additional dynamo action beyond pure compression. 
We note that we find stronger magnetic field amplifica- 
tion for the sub- than for the super-isothermal situation. 
However, the interpretation of this result is difficult, be- 
cause for T < 1 the collapse is filamentary and so the 
normalization by a factor p 2 / 3 does not compensate com- 
pletely for the contribution from flux-freezing. 

4. DISCUSSION AND CONCLUSIONS 

We performed a set of four high-resolution numerical 
simulations, where we studied gravitational collapse and 
magnetic field amplification for four different thermody- 
namic regimes. We find that the dynamics of gravita- 
tional collapse and consequently the star formation pro- 
cess depends significantly on the thermal evolution of 
the gas (e.g. Larson 1985[ ). If the gas is able to heat 
up during compression, i.e. for an effective T > 1, we 
find the formation of a central virialized core, leading 
to the production of turbulent motions and highly tan- 



gled magnetic field stru ctures (see also Sur et al. 2010 
Federrath et al. 2011b). Under sub-isothermal condi- 
tions (F < 1), the formation of such cores is suppressed. 
The collapse proceeds much faster and the flow develops 
strong shocks. These create a network of intersecting 
sheets and extended filaments. In this case, vortices still 
form in the vicinity of the c urved shocks and at the inter- 
section of multiple shocks ( |Sun fc Takayam a 2003), but 
are absent in the infalling region. 

We find that magnetic fields are amplified in the super- 
isothermal as well as in the sub-isothermal region. The 
local magnetic field strength is increased both by com- 
pression, which is more efficient within the shocks, and by 
the presence of turbulent motions. The turbulent mag- 
netic field amplification can be described in terms of the 
Kazantsev theory, and has been shown to depend on the 
square root of the Reynolds number for Kolmogorov tur- 
bulence and on the cubi c root of the Reynolds number 
for Burgers turbulence ([Kazantsev 1968[ Subramanian 
1998 Schober et al. 2012] ^ The transsonic simulations 
presented here lie in between these two extreme cases. 
Because of the large Reynolds number of astrophysical 
gases, the inferred amplification timescales are typical of 
order of, and under most condit ions significantly shorter 



than, the local free-fall time (Federrath et al. 2011b 



Schober et al.||2012). We therefore expect the presence 



ol magnetic fields close to the saturation level for a wide 
range of T (even for T < 1) quite independent of environ- 
mental conditions. In particular, we argue that magnetic 
fields must not be neglected when studying stellar birth 
in the high-redshift Universe. 

We expect the equation of state to play a major 
role also under different conditions, for instance in the 
observe d filamentary struc t ure of nearby star- formin g 



clouds (Andre et al. 2010 Arzoumanian et al. 2011). 



We note that the interstellar medium (ISM) in the so- 
lar neighborhood at number densities n in the range be- 
tween 10 and ~ 10 4 cm -3 lies in a strong cooling regime 



with r « 0.7 ( jLarson|i 985), where cooling i s provided by 
the lin e emission mostly of C + and CO ( |Glover fc Clark| 
|2012a| ), before coupling to dust results in F 



densities above n ~ 10 cm 



1.1 for 

It is well established that 



the supersonic turbulence ubiquitiously observed in the 
ISM by i tself will lead to a highly fi lamentary morphol- 
ogy (e.g. Mac Low fc Klessen||2004[ ). Here we want to 
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Fig. 3. — Velocity structure for the runs with F = 0.7,0.9 (from top to bottom) at the time shown in the first column. The columns 
show density, tff divv and % (curl v)^, respectively. The arrows are velocity vectors. The velocity mentioned in the images corresponds to 
an arrow of the length given in the legend. The green circle marks the Jeans volume. 
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Fig. 4. — Velocity structure for the runs with F = 1.1, 1.2 (from top to bottom) at the time shown in the first column. The columns 
show density, % divv and £ff (curl v)z, respectively. The arrows are velocity vectors. The velocity mentioned in the images corresponds to 
an arrow of the length given in the legend. The green circle marks the Jeans volume. 
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Fig. 5. — Evolution of the turbulent velocity dispersion, effective 
magnetic field amplification, ratio between the small-scale mag- 
netic field and the turbulent kinetic energy, and the turbulent Mach 
number within the Jeans volume for the collapse simulations with 
r = 0.7, 0.9, 1.1,1.2. 

draw attention to the fact that also the thermodynamic 
response of the gas plays an important role in shaping 
the p hysical chara cteristics of the observed filaments (see 
also [Larson 2005). We argue that the cooling processes 
associated with the formation of molecular clouds are a 



key to understand the observed morphological structure 
of nearby star- forming regions. 
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